Operator Product Expansion for Exclusive Decays: B + — > D+e + e and B + — > D* 
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The decays B + — > D + d e + e~ and B + — ► D* d e + e~ proceed through a weak and an electromagnetic 
interaction. This is a typical "long distance" process, usually difficult to compute systematically. We 
propose that over a large fraction of phase space a combination of an operator product and heavy 
quark expansions effectively turns this process into one in which the weak and electromagnetic 
interactions occur through a local operator. Moreover, we use heavy quark spin symmetry to relate 
all the local operators that appear in leading order of the operator expansion to two basic ones. We 



use this operator expansion to estimate the decay rates for B 



D 



(*)+„+, 



£>-mesons, with all their different decay channels, have 
become one of the prime objects for checking the stan- 
dard model and measuring its parameters. This has led 
to the development and ongoing planning of facilities and 
experiments dedicated to their study. In the near future 
hadronic facilities will produce more than 10 11 B- 
mesons per year, allowing studies of very rare processes. 
Only processes which afford a trustworthy quantitative 
theoretical description can be used to verify the stan- 
dard model. There has been some interest recently in the 
decays B+ -► L>*+ 7 @ and B+ -> D d *l + e+e- g. How- 
ever, it was already realized in [^]||] that there were tech- 
nical problems with the calculations. We demonstrate 
that over a large fraction of phase space a combination 
of an operator product and heavy quark expansions ren- 
ders B + — > D d e + e~ computable. Our methods are 
fairly general and can be applied to a variety of other 
processes. 

The Operator Product Expansion (OPE) is commonly 
used in the calculation of inclusive decay rates. One uses 
the optical theorem and the performs an OPE on forward 
scattering amplitudes. We will show that for a class of ra- 
diative decays one may use the OPE to compute exclusive 
decay amplitudes. There is a simple physical motivation 
for the use of the OPE in the decay amplitude of, for 
example, B + — » Z?(*)+e + e~ . Consider the hadronic part 
of the amplitude to lowest order in electromagnetic and 
weak interactions, 



(1) 



Here O = 67 l/ (l — 75)1x07^(1 — 75)^ is a four quark opera- 
tor responsible for the weak transition, j£ m is the electro- 



magnetic current and T stands for time ordering of these 
operators. The energy denominators from intermediate 
states of energy E are Mb — E; since the B mass, Mb, is 
much larger than that of the D meson and its excitations, 
many intermediate states contribute significantly to the 
amplitude. Since the available energy is much larger (in 
units of Aqcd) than the energy spacings between inter- 
mediate states, the time ordered product in (^) should 
be well approximated by an expansion in local operators 
(OPE). 

This OPE is not a short distance expansion since 
the momentum transfer is not in the Euclidean domain. 
Hence its validity relies on quark-hadron duality. This is 
exactly analogous to the the use of an OPE in the com- 
putation of heavy hadron lifetimes. For large enough Mb 
violations to duality can be described by (uncomputable) 
powered suppressed operators in the OPE Q. Therefore, 
we can trust at least the leading term in our computation 
of the amplitude. 

To use an OPE in Eq. ([[]) we take the heavy quark 
limit for the b and c quarks. Rather than the matrix el- 
ement with external physical mesons of Eq. (|l|), we first 
consider a Green function with two external quarks and 
two external anti-quarks. For the momenta of the heavy 
quarks we write rribV + k\, and m c v' + k c , while for the 
light quarks we take k u and kd- The residual momenta 
are small, ki <C ?ri(, jC , provided w = v ■ v' remains of order 
unity. Just as in the case of semileptonic inclusive de- 
cays H the combined heavy-quark and operator product 
expansions give an expansion in powers of ki/m^c- 

To see how this works explicitly, consider the first 
Feynman diagram of Fig. |l|. To lowest order in fej one 
has 



m b j, + m c 



7"(1 -Ts) <8 7i/(l -75), 



(2) 
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where Q c = 2/3 is the charge of the c quark. This corre- 
sponds to a local operator which depends on the heavy 
masses and the velocity v, which is not a kinematic vari- 
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u 






FIG. 1. Feynman diagrams representing lowest order con- 
tributions to the Green function. The filled square represents 
the four quark operator O and the cross represents the elec- 
tromagnetic current cf. Eq. (|l|). 



able but a parameter in the heavy quark effective theory. 
Corrections appear as higher dimension operators sup- 
pressed by powers of the large mass. For example, the 
leading correction is of the form 



7 (1-75) ® 7i/(l -75). 

Of course, in the operator language, k c becomes —id act- 
ing on the c-quark field. 

Consider next the second Feynman diagram of Fig. [I]. 
Neglecting the light quark masses we have 

-Q^(l-7 5 ) r W " Wc? !; 2 7^>(l-7 5 ). (3) 
{m b v — m c v ) z 

This term differs in an important way from (||). The 
large denominator is not uniformly large over the whole 
kinematic range. In fact the denominator is just the 
square of the leading part of the momentum out of the 
electromagnetic current, q = m b v — m c v' + ^ki. For 
the decay B + — > £>(*' + e + e~, the kinematic range is 
< q 2 < (Mb — Mj)) 2 . The approximation is valid 
provided Aqcd *C m Cib , e.g., the corrections are or- 
der Aqcd/toc,6- There are also corrections of order 
h-QCD m b.c/ q 2 ■ So our results are limited to the region 
were q 2 scales like m 2 b . The region were q 2 does not 
scale like m 2 b is parametrically small, so the arguments 
we present are theoretically sound. We emphasize that 



our method cannot be applied to B + — > D + j, for which 
q 2 = 0, identically. 

The third and fourth diagrams in Fig. |l| are similarly 
computed. In an obvious notation, these are 



,(1 - 75 ) ® 7"-^ ^7^(1 " 75) 



and 



Qbl, 



Q^(l-75)^7-(l-7 5 ) r W " m ^ 2 7^ 

{m b v — m c v y 



(4) 



(5) 



Finally, for the calculation of the rate B + — > D^*' + e + e" 
one must compute the matrix elements of these four local 
operators. 

The matrix elements of all of these local operators can 
be expressed, by use of heavy quark spin symmetry, in 
terms of only two invariant functions. Spin symmetry 
is best exploited using the Wigner-Eckart theorem. We 
represent the meson by 



£>75 



and the anti-meson by 



h® = [b:y - S75] 



1 - j> 



We will also need the field conjugate, 



i+f 



(6) 



(J) 



(8) 



In terms of these we find 
(H^\h^r c dh^T b u\H^) = 

^-Tr(B^T c )Tr(H^T b ) + ^Tr(H^T c H^T b ), 

(9) 

where T b ^ c are arbitrary 4x4 matrices, h v , is the field 

that creates a heavy quark with velocity v' and hv an- 
nihilates a heavy anti-quark with velocity v. Invariance 
under spin symmetry of the heavy quarks readily implies 

{H® \h$T c d h^T b u\H^) « H$>T e ® H®T b . (10) 

Invariance under rotations implies that the remaining 
four indices must be contracted. 

The functions (3 and 7 can be determined in simula- 
tions of QCD on the lattice. In the calculations below 
they are estimated using a vacuum insertion approxima- 
tion, which gives (3(w) = /bId^/MeMd and "f(w) = 0. 
A better estimate can be obtained by using isospin and 
heavy flavor symmetry to relate (|^) at w = 1 to 



2 



(fl|67"(l - 75)^(1 - j 5 )d\B) = lB B f B M%, (11) 
(B\b(l - l5 )db(l - j 5 )d\B) = -§B s f B M 2 . (12) 

Isospin relates these matrix elements to the matrix el- 
ement in (^) symmetrized in u «-> d. This symmetrized 
operator can be related || to (^|) if we make the following 
assumption: 



H^\Jt)T c T A d Wr b T A u\H^) = 0. 



(13) 



Here T A is a generator of color-S'C/(3). This assump- 
tion is weaker than vacuum insertion because it holds 
true even if we insert a complete set of physical states, 
rather than only the vacuum, between the currents. We 
obtain 13(1) = \(B B + bB s ) f B f D y/M B M D and 7 (1) = 
~l{B B ~B s )f B f D ^M B M D . A recent calculation us- 
ing quenched lattice QCD gives B B ss B$ ~ 0.8, with a 
strange light quark and a c and b heavy quark for B B and 
Bs, respectively. Therefore our calculations below, which 
uses the vacuum insertion values (3(1) — f B foVM B MD 
and 7(1) = can be trivially modified to account 
for these lattice results, (3(1) « 0.8/b/d\/Mb^d and 
7(1) ~ 0. We do not use these lattice results below 
because B B and Bs are computed for different heavy 
quarks and different renormalizations, and the result for 
B B differs by 20% from other lattice studies 

The effective Hamiltonian for the weak transition is 

H' eS = ^ V ub V c * d (c(fi/M w )0 + c 8 ( M /M w )0 8 ) , (14) 

where C 8 = 67^(1 - j 5 )T A ucry v (l - ~f 5 )T A d. The de- 
pendence on the renormalization point \i of the short 
distance coefficients c and c 8 cancels the /x-dependence 
of operators, so matrix elements of the effective Hamil- 
tonian are /x-independent. Resuming the leading logs, 



c(fi ) = |x 2 + |x 1 and cgQ^o) 
x = (a s (pL Q ) /a s (M w )) 6/23 . 
Defining 



where 



h (*)M = (jD (*)+| / # x e *x T(j£ m (x)H e g(0))\B + ), (15) 



the decay rate for B + — » D^ + e + e is given in terms of 



q 2 and t = (po - 
dT 



Pe+, 



{pb -Pe-r by 



dq 2 dt 



2 8 7T 3 M| 



2 V 



(16) 



where = u(p e -)j^v(p e +) is the leptons' electromag- 
netic current. A sum over final state lepton helicities, 
and polarizations in the D* case, is implicit. Using the 
OPE and spin symmetry we obtain 



(Awm b — 3m c )u^ — {Zm b — 2wm c )v'^ 



(m b v — m c v') 2 
m b v lfl 



m c v h 



IDT, 



rat 



(17) 



and 



m b (3e^ - 4u • ev^) + m c (v ■ ev'^ - 3we M ) 



(m b v — m c v') 2 



(m b v — m c v') 



A2 



(18) 



mfce M + m c (v ■ ev lfl — we^) + im c e^ a,3 ~> 'e a v'gVj 



+ 



Here k = G F /V^ V ub V* q [c(3 + c 8 (3 8 ] and we have defined 
the matrix element of 0$ as (3s, in analogy to Eq. (||). 
These expressions are our central results, demonstrating 
that the decay rates for B + — > D^ + e + e~ can be ex- 
pressed in terms of the matrix elements (3 and (3g. Below 
we make an educated guess of these matrix elements, 
but for reliable results they should be determined from 
first principles, say, by monte carlo simulations of lattice 
QCD. 

There are short distance QCD corrections to the coef- 
ficients (0)-(||) in the OPE. In leading-log order they are 
determined by the renormalization group. The matrix 
elements of the operators are all related by spin sym- 
metry to the reduced matrix elements (3, (3%, 7 and 7s- 
However, (3 and (3% do not mix into 7 and 7s. Moreover, 
since we use 7 = 73 = 0, we need only the 2x2 matrix 
of anomalous dimensions of (3 and (3g . A straightforward 
one loop computation gives the anomalous dimension for 
(3 and (3s as follows: 



7 : 



8 4(wr(w) - 1) 

4^ V fKH - 1) l(7-wr(w)) J ' 



(19) 



where r(w) = . \ ln(u> + \Jw 2 — 1). Note that the 
diagonal entry for (3 is precisely twice the anomalous di- 
mension of the heavy- light current, as if the operator 
was factorized. This is in accordance with the results of 
calculations of 7 of the operator for B — B mixing ||. 
Moreover, at w = 1 the matrix simplifies, 7 = =2*1. 

Denote by c and c$ the coefficients of (3 and respec- 
tively. They satisfy a renormalization group equation, 

* 5 = _ 7 t £ . (20) 
u/i 

Here "T" denotes transpose of a matrix and c is a column 
vector, c T — (c, eg). In leading-log order the solutions are 

c( M ) = z^[\(2zJ + z-t)c(jM>) + l(zt- z-t) c 8 (mo)] 



f'8 



(/i) = z+[(zt ~ z"«) c(/xo) + i + 2z~t) c 8 (/xo)] 



where z = a s (/x)/a s (/j ), 
13 — wr(w) 



TP 



wr(w) — 1 



Sb 



(21) 



b a is the coefficient of the one loop beta function in QCD, 
bo = 11 — with rif = 3 light flavors in our case. 
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For our numerical estimates below we match the coef- 
ficients c and C8 to c and c 8 at the intermediate scale 
Ho = y?mi,m c . Numerically, with a s (/j,o) / a s (Mw) ~ 2.24 
one has c « 1.0 and c 8 » —0.7. We will use fi = 1 GeV, 
with the implicit understanding that the matrix elements 
(3{w) and (3%{w) are computed at that value of the renor- 
malization point. 

It is now a trivial exercise to compute the differential 
decay rate. We integrate the rate in Eq. (|l6|) over the 
variable t and obtain, for B + — > D\ 



dT 



dq 2 288ttM| 



|K b V cs | 2 (c/3 + c 8 /3 8 ) 2 ^(g). (22) 



Here T(q) is a dimensionless function of q = ^/ q 2 jm 
and 771 = MoJMb- For £? + — > Z?+e + e~ it is given by 

4(2 - m 2 ) 2 [(l - (<? + m) 2 )(l - (<? - m) 2 )] 3 / 2 



3<7 4 7n(l — m 2 ) 5 



of heavy mesons. At zero hadronic recoil the matrix ele- 
ment of the leading operator in the OPE is related by 
heavy quark symmetries and octet suppression to the 
matrix element for BB mixing. Although the rates we 
compute are too small to be observable at B-factories, 
they may be accessible to experiments in hadronic col- 
liders. We have demonstrated the method by calculat- 
ing the rate for B + — > £r ' d e + e~ to leading order in 
the operator expansion and to leading- log order in QCD. 
Systematic corrections to both of these expansions could 
and should be computed. Irreducible errors are expected 
from the application of the operator product expansion 
in the time-like regime, which is analogous to the as- 
sumption of local quark-hadron duality in calculations of 
heavy hadron lifetimes. 

Acknowledgments This work is supported by the De- 
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97ER40546. 



(23) 



and for B 



D*+e+e- by 



4 x /(l-(g + rn) 2 )(l-(g-rn) 2 7 
3 g 6 m(l — m 2 ) 2 



(-36m 8 + rh 2 q 8 



«fm 4 



-17m 8 q 2 



30g 4 m 4 



9m 10 - 
38g 2 7h 6 



t- 8m 6 q 4 
4q 6 m 2 + 4<f 



+4<f - 8q 4 - 36m 4 + 9m 2 - 4<f m 2 
+30<fm 2 - 21g 2 m 4 + 54m 6 ). 



(24) 



Note that we have not distinguished between heavy quark 
and meson masses. The distinction enters at order I/tuq 
in the heavy mass expansion, and we have not considered 
such corrections in this work. 

To estimate the branching fraction numerically, we 
use the vacuum insertion approximation, (3{w) — 
z~ A/bo .fBfDVM B M D and fo(w) = 0, and integrate the 
rate from q 2 = 1 GeV up to the kinematic limit q 2 = 
(Mb — Mo s ) 2 ■ The lower limit is an estimate of how low 
q 2 may be before our operator expansion breaks down. 
We find 
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-11 


(27) 


[7 


Br(J3 + - 


-> D+e + e- 


~)\q 2 >l 


GeV 


= 1.4 x 10" 


-11 


(28) 



rgi 

where we have used IKsKal = 0.004, i KbKrfl = 8.8 x 
1(T 4 , f B = 170 MeV and f D = f B ^M B /M D . It is 
important to observe that the portion of phase space q 2 > [9] 
1 GeV is expected to give a small fraction of the total rate 
since the the pole at q 2 = dramatically amplifies the 
rate for small q 2 . 

In summary, we have applied the operator product and 
heavy quark expansions to the exclusive decay amplitude 
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